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ABSTRACT 

Diffuse interstellar clouds show large abundances of HJ which can be maintained only by a high 
ionization rate of H2. Cosmic rays are the dominant ionization mechanism in this environment, so 
the large ionization rate implies a high cosmic-ray flux, and a large amount of energy residing in 
cosmic rays. In this paper we find that the standard propagated cosmic-ray spectrum predicts an 
ionization rate much lower than that inferred from . Low-energy 10 MeV) cosmic rays are 
the most efficient at ionizing hydrogen, but cannot be directly detected; consequently, an otherwise 
unobservable enhancement of the low-energy cosmic-ray flux offers a plausible explanation for the 
it^ results. Beyond ionization, cosmic rays also interact with the interstellar medium by spalling 
atomic nuclei and exciting atomic nuclear states. These processes produce the light elements Li, Be, 
and B, as well as gamma-ray lines. To test the consequences of an enhanced low-energy cosmic- 
ray flux, we adopt two physically-motivated cosmic-ray spectra which by construction reproduce the 
ionization rate inferred in diffuse clouds, and investigate the implications of these spectra on dense 
cloud ionization rates, light element abundances, gamma-ray fluxes, and energetics. One spectrum 
proposed here provides an explanation for the high ionization rate seen in diffuse clouds while still 
appearing to be broadly consistent with other observables, but the shape of this spectrum suggests 
that supernovae remnants may not be the predominant accelerators of low-energy cosmic rays. 
Subject headings: astrochemistry - cosmic rays - ISM: clouds - ISM: molecules 



1. INTRODUCTION 

Several recent observations of Hj^ in the diffuse 
interstellar medium (ISM) indicate an average cosmic- 
ray ionization rat e of molecula r hyd ro gen of about 
4x 10-16 s-i (|McCafl et al.l [200I llndriolo et aP 
l2007t ). This value is about 1 order of magnitude larger 
than was previously inferred us ing other m olecular 
trace rs such as HP and OH (jO'Donnell fc Watson 
1971: IMadt fc Dalgarnd MtA iBlack etal.l HOTSl : 
Hartauist et al.l Il978al lbl: iFederman et al.l |1996[ ). How- 



ever, several models have also required ionization rates 
on the order of lO'^^ s"^ Ivan Dishoeck & Black '1986'; 
iLisztl 120031 : iLe Petit et all [2004 : Shaw et al. .2006.. ,2008) 
in order to reproduce the observed abundances of 
various atomic and molecular species. This agreement, 
coupled with the simplicity behind the chemistry of 
Hj^, leads us to conclude that the newer measurements 
are most likely correct. In this paper, we explore the 
implications that a high ionization rate has for cosmic 
rays and related observables. 

Aside from observational inferences, the cosmic-ray 
ionization rate can also be calculated theoretically us- 
ing an ionization cross section and cosmic-ray energy 
spectrum. While the io nization cross section for h ydro- 
gen is wefl determined ()Bethelll933t IInokuti|[l97ll ). the 
cosmic-ray spectrum below about 1 GeV is unknown. 
This is because low energy cosmic rays are deflected 
from the inner solar system by the magnetic field cou- 
pled to the solar wind (an effect called modulation) 
and so the flux at these energies cannot be directly ob- 
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served. This theoretical calculation of t he ionization rate 
has been perfo rmed several times (e.g., iHa yakawa et al.| 
1£61; Spitzcr fc Tomaskolll968t iNath fc B iermann 199$ 
Webber! 1199 8). with each study choosing a different low 
energy c osmic-ray spec trum for various reasons. Most 
recently, IWebbeij (|1998l ) predicted an ionization rate of 
(3 — 4) X IQ-i^ using interstellar proton, heavy nuclei, 
and electron cosmic-ray spectra. The proton and heavy 
nuclei spectra were found by attempting to remove the 
effects of solar modulation from Pioneer and Voyager ob- 
servations, while the electron spectrum was derived from 
radio and low-energy gamma-ray measurements. Even 
having accounted for all of these components, this result 
falls about 1 order of magnitude short of the inference 
based on , suggesting that the de- modulated solar 
system spectrum may not be the same as the interstel- 
lar spectrum in diffuse clouds, and/or that the IWebbed 
(1998) extrapolation to low energies underestimates the 
true interstellar value. 

Together, all of the above studies have shown that 
the ionization of interstellar hydrogen is a powerful 
observable for probing cosmic-ray interactions with the 
environments through which they propagate. Beyond 
ionization though, cosmic rays will interact with the ISM 
in other ways which lead to additional and complemen- 
tary observables. Namely, inelastic collisions between 
cosmic-rays and interstellar nuclei inevitably: (i) create 
light element isotopes ^Li, ''Li, ^Be, ^^B, and ^^B when 
cosmic rays spall C, N and O nuclei ([Reeves 
Meneguzzi ct aT[ |1971[ ). and (ii) excite nuclear states 
such as i^C* and ^^O , the decay of which produce 
gamma-ray lines, most p rominently at 4.44 MeV an d 
6.13 MeV, respectively (i Meneguzzi fc Reeves! Il975al) . 
Similar to the theoretical calculation of the ionization 
rate, a cosmic-ray spectrum and relevant cross sections 
can be used to determine the production rates of light 
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2. THE IONIZATION RATE INFERRED FROM H+ 

The chemistry behind H;^ in the diffuse ISM is rather 
simple. Its formation and destruction are given by the 
reactions: 



CR + Ha ^ CR 



H+ + e- 



H 



H, 



H2 + H or H + H + H. 
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Fig. 1. — Example cosmic-ray spectra used in t h e lite rature 
over the past few decades. 1 - 'Herbst & Cuppen' f2006'); 2 - 
[Spitzcr & Tomasko (1968); 3 - Gloccklcr & Jokioii (1969) (via fit- 
ting function of j^ncUcr ct al. (2003)); 4 - Id & Axford (1985) (via 
fitt ing function of Vallc ct al. (2002)); 5 - Havakawa ct al. (1961); 
6 - INath &: Biermann (,1994 ). The dotted line is the propagated 
leaky box spectrum used in this paper, also shown in Figure |2] 
Note the agreement above and discrepancy below 1 GeV. These 
spectra were selected to be illustrative of choices in the literature 
used for diffe r ent a pplications. Of these, it is worth noting that the 
Up &: AxfordI l|1985l ) spectrum is the result of a calculation specifi- 
cally designed to recover the (demodulated) propagated interstel- 
lar cosmic-ray spectrum. The shaded region shows the range of 
uncertainty in the demodulated proton spectrum as described by 
IMoril (119971 ). Crosses arc proton data from the top of the Earth's 
atmosphere (AMS Collaboration ct al. 2002 ) and clearly show the 
effects of modulation below ~ 1 GeV. 



H2 is first ionized, after which the ion quickly re- 
acts with another H2 molecule to form . In diffuse 
(and dense) clouds it is assumed that this ionization is 
due almost entirely to cosmic rays, as the flux of pho- 
tons with E > 13.6 eV will be quickly attenuated by 
atomic hydrogen in the outer regions of the cloud. The 
first step is many orders of magnitude slower, so it acts 
as the overall rate limiting step. The primary channel 
for destroying H^j" in diffuse clouds is recombination with 
an electron. H;^ is destroyed by reaction ([3]) on a time 
scale of about 100 years, much shorter than the~ 10^ 
lifetime of diffuse clouds ( Wagenblast & Hartq uistlllQSS 
so the steady-state approximation is valid and the forma- 
tion and de struction rates can b e equated. This assump- 
tion yields (jGeballe et al.lll999[ ) 



C2n(H2) - fce«(H+)n(e), 



(4) 



where C2 is the ionization rate of H2, fee is the H;^-electron 
recombination rate constant, and ?i(X)'s are number den- 
sities. Spectroscopic observations of transitions from 
the two lowest rotational levels of H;^, the only levels 
populated at the low temperatures of diffuse interstellar 
clouds, provide the column density. The cloud path 
length is then found using the observed hydrogen column 

^ . density and inferred hydrogen number density. Divid- 

ment / e .g., iMeneeuzzi et al ] 1 197 1'; 'Meneguzzi & Reeves' ing the column by the path length gives n{'R'^), and 

leaves three variables in the steady state equation: ke, C2j 
and n(e)/n(H2). However, previous work has shown that 



elements a nd gamma-ray lines. Both the light ele- 
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(e.g., [M eneguzzi & Reeves 1975a; Ramatv et al. 1979 
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2iiiTir calculations have been performed multiple times, 
again with each study choosing a different low energy 
cosmic-ray spectrum. 

Some of the spectra that have been used for these cal- 
culations are shown in Figure [H While several more 
cosmic-ray spectra have been used, many share func- 
tional forms with those plotted and so have been omitted 
for the sake of clarity. Note that all of the spectra agree 
with demodulated data (shaded region) above a few hun- 
dred MeV and raw data (crosses) above a few GeV, but 
that they can differ by about 4 orders of magnitude at 1 
MeV. Figure [1] is shown primarily to illustrate our poor 
understanding of the low-energy portion of the cosmic- 
ray spectrum. 

In this study, we calculate the cosmic-ray ionization 
rate using several low energy cosmic-ray spectra in an 
attempt to reproduce the value inferred from obser- 
vations. For the spectra that successfully predict ioniza- 
tion rates close to the inferred value of 4 x 10""'^^ s~^, 
we further investigate the implications that they have on 
dense cloud ionization rates, light clement abundances, 
gamma-ray fluxes, and energetics arguments. 



the H j"-electron recombination rate (jMcC all et al."200 



20041) and the electron-to-hydrogen ratio (CardcUi et a! 
1996( 1 are relatively well constrained, leaving ^2 as the 



only free parameter. Starting from eq. ^ and using var- 
ious o ther relationships and assumptions. Ilndriolo et aTl 
(|2007f l derived an equation for the cosmic-ray ionization 
rate that depends on observables. This equation was 
then used to infer the ionization rate toward several dif- 
fuse cloud sight lines. From all of the sight lines with 
Hj^ detections, the average cosmic-ray ionization rate of 
molecular hydrogen was found to be about 4 x 10""'^^ s~^ 
with a maximum unc ertainty of about a fa ctor of three 
either way (see §4.2 of Ilndriolo etHl (|2007D for a discus- 
sion of the calculations and uncertainties) . 

3. IONIZATION ENERGETICS: A MODEL-INDEPENDENT 
LOWER BOUND 

Assuming that the ionization rate above is uniform 
throughout the diffuse Galactic ISM, it is relatively sim- 
ple to estimate the total. Galaxy-wide amount of power 
necessary to produce such a high value. While this as- 
sumption of uniformity is not strictly valid (there are 
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and clouds in the 



diffuse clouds with (2 

Galactic center with C2 >, 10'^^ s~ \ (lOka et al.l l2005t 
lYusef-Zadeh et al.l l2007t iGoto et al.l l2008f )). the small 
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sample size of sight lines does not allow for the determi- 
nation of a meaningful relationship between position and 
ionization rate. Despite these fluctuations, if all atomic 
hydrogen experiences the same ionization rate on aver- 
age, then the Galactic luminosity in ionizing cosmic rays, 

.ionize? 

is given by 



T /■ 'TTF f -^^H, diffuse 
-t^CR, ionize = C,H Aii 

V ma 



(5) 



where AE is the average energy lost by cosmic rays per 
ionization event. The number of hydrogen atoms in dif- 
fuse clouds is the ratio of the mass of all atomic hydrogen 
in diffuse clouds in the Galaxy, Mh, to the mass of a hy- 
drogen atom, mn- The ionization rate of atomic hydro- 
gen, (^H, is related to the ionization rate of molecular hy- 
drogen (the observable pro bed by H^) by 2.3Ch = 1-5C2 
(|Glassgold fc Langeilll97^ . The coefficients here are fur- 
ther explained in §5. 

Given the ionization rates from the previous section, 
we may place a model-independent lower limit on the 
ionizing cosmic-ray luminosity as follows. Each ioniza- 
tion event requires a cosmic-ray energy input AE > 13.6 
eV, the ionization potential of atomic hydrogen. On av- 
erage cosmic rays will lose more than 13.6 eV in the ion- 
ization process though, so by setting AE — 13.6 eV in 
eq. (O we calculate a hard lower limit on the power re- 
quirement. The gas mass relevant to eq. ([5]) is that of 
all neutral hydrogen in Galactic diffuse clouds, which 
we take to be MndifFuse = = 5 x lO^Mp^ (the average 
of IHenderson et all (119821) . ISodroski et all 11994!), and 
iMisiriotis et al.l (|2006l )). This value results in a lower 
limit to the cosmic-ray luminosity of 



LcRaonize > 0.11 X 10^^ erg (100 yr)-i 



^'^H, diffuse 
5 X lO^M© 



Note that this cosmic-ray "energy demand" is in addi- 
tion to the requirements found based on co smic-ray en- 
ergy l ost as the particles escape the Galaxy. iFields et al.l 
(|2001[ ) estimated the sum of both contributions, and 
found LcR,tot — 0.5 x 10^^ erg (100 yr)~^ which is con- 
sistent with eq. ([6|) but also implies that ionization repre- 
sents a significant part of the cosmic-ray energy budget. 

However, eq. ^ is only the lower limit to the amount 
of cosmic-ray energy that goes into ionization. We can 
get an actual estimate on the luminosity of ionizing 
cosmic rays by accounting for molecular hydrogen and 
by using a more precise v alue of AE. According to 
ICravens fc Dalgarnol (|1978D the average energy lost dur- 
ing an ionization event is about 30 eV, which by itself 
increases LcR.ionize to 0.24 x 10^^ erg (100 yr)~^. The 
inclusion of H2 is more complicated. The mass of H2 is 
about the same as that of H in the Galaxy, but m ost H2 
resides in dense molecular clouds (|Brinksl 1199 01) which 
do not experienc e the same cosm ic-ray ionization rate as 
the diffuse ISM (|Dalgarnol [2006t ) . Assuming half of aU 
Galactic H2 experiences the ionization rate used above, 
-C'CR, ionize ~ 0.33 X 10 ^^ erg (100 yr)~^, a large fraction 
of the result found bv lFields et"an (|200lD . 

As it is currently believed that Galactic cosmic 
rays are accelerated in supernova remnants (SNR), 
these results have implications for the efficiency with 
which supernova mechanical energy is transferred to 



particle acceleration. If a typic al supernova releases 
lOf^erg of mechanical energy (e.g.. I Arnettll 19871 : 1 WooslevI 
|1988( ) and 3 ± 2 s upernovae (SNe) occur each cen- 
tury in the Galaxv (|van den Bergh fc TammannI 119911 : 
iDragicevich et al.l [19990, then at least 4% of the energy 
released in SNe must accelerate the cosmic rays which 
ionize hydrogen in the ISM. This efficiency climbs to 12% 
if we take the more realistic estimate instead of the lower 
limit. However, uncertainties in the supernova rate, su- 
pernova energy, and mass of hydrogen in the Galaxy lead 
to a total uncertainty of about a factor of 5 either way for 
this value. It is important to note though that this cal- 
culation depends only on the cosmic-ray ionization rate, 
and not on an adopted form of the cosmic-ray spectrum. 
In contrast, calculating the ionization rate is highly de- 
pendent on the cosmic-ray spectrum, to which we now 
turn. 

4. POSSIBLE SPECTRA OF LOW-ENERGY COSMIC-RAYS 

Given the well-understood physics of the passage of 
energetic particles through matter, the ionization rate 
completely reflects the spectrum of cosmic rays. In par- 
ticular, the ionization cross section (below, eq. I12p grows 
towards low energies as tJion ~ v^'^ ^ which means 
that the lowest-energy particles have the strongest effect 
on ionization. Given our lack of direct observational con- 
straints on cosmic rays at low energies, we will examine 
the ionization arising from various possible low-energy 
behaviors which are physically motivated and/or have 
been suggested in the literature. Here we summarize in 
a somewhat pedagogical way some of the main features 
of the current understanding of cosmic-ray acceleration 
and propagation. 

The cosmic-ray spectrum with the strongest physical 
motivation (in our view) takes supernova explosions to 
be the engines of Galactic cosmic-ray acceleration. That 
is, supernovae remnants provide the sites for diffusive 
shock acceleration and thus act as cosmic-ray sources. 
At these sources, diffusive shock acceleration creates par- 
ticles with spectra which are close to simple power-laws 
in (relativistic) momentum p. Specifically, consider the 
"test-particle" limit when particle acceleration has a neg- 
ligible effect on the shock energy and structure. In this 
limit, the cosmic-ray production rate, q, per unit volume 
and time and per unit interval in relativistic momen- 
tum has famously been analytically shown to be (e. g., 
lKrvmskiil[T977l : lBeIl[l97at iBlandford & OstrikedlT97^ 



Qshock £ 



dV dt dp 



oc p 



(7) 



i.e., a power-law in momentum. Here the momentum in- 
dex in the case of strong shocks is x — 2-t-4/A^^, where 
the shock Mach number is A4 = Wshock/csound,ism ^ 
1. The upshot is that for strong shocks (large A^) as 
one would find in supernova remnants, the acceleration 
power law is just slightly steeper than the flattest (i.e., 
largest at high-energy) limiting power-law spectrum al- 
lowed by energy conservation: qum cx p~^. Going be- 
yond the test-particle limit requires nonlinear treatment 
of the feedback of cosmic ray energy and pressure on the 
shock structure and evolution; the study of this nonlin- 
ear shock acceleration remains a vital field, but several 



group s (e.g.. iKang fc Joned 119951: iBerezhko fc^Uison 
Il999t iBlasi I2002D find that the accelerated particles have 
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a spectrum which is roughly similar to the test-particle 
result, but which shows some concavity in momentum 
space, i.e., the effective spectral index x — dluq/ d\n.p 
does show a modulation around the constant test-particle 
value. Intriguingly, there seems to be agreement on the 
qualitative result that the low-energy flux will be higher 
than for the test-particle predictions. Unfortunately, 
the quantitative results remain at present rather model- 
dependent. For the purposes of our analysis, we will sim- 
ply adopt test-particle power-law acceleration spectra as 
in eq. ([7]). Our results can then be viewed as testing the 
validity of the test-particle approximation at low ener- 
gies. 

Once produced at acceleration sites, cosmic rays prop- 
agate away, and eventually are removed either by es- 
cape from the Galaxy or by stopping in the ISM due 
to energy losses (predominantly by energy transfer to 
the ISM, either ionization or excitation of atoms or 
molecules). Propagation alters the spectra of cosmic 
rays from those at the sources. Theoretical treatments of 
these effects typically make the simplifying assumption 
of a steady state balance between production and losses. 
The resulting "propagated" spectrum should represent 
the flux as seen by an average region of the interstellar 
medi um, far from cos mic-ray sources (elegantly reviewed 
in .Strong et al.ll2007[ ). 

A full calculation of cosmic-ray propagation at mini- 
mum involves the particle "flows" in energy space; the 
simplest such treatment is the classic "leaky-box" model 
which treats the Galaxy as a medium with sources dis- 
tributed homogeneously. More sophisticated models ac- 
count for the inhomogeneous Galaxy and effects such as 
diffusion and re-acceleration. In ge neral, when models in- 
clude the low-energy regim e (e.g., iLerche fc Schlickeised 
Il982t IShibata et all l2006f) . they find that when ini- 
tially accelerated or "injected" spectra are power-laws 
in momentum, the resulting propagated spectra are very 
nearly also power laws, with fairly abrupt changes of 
spectral indices at characteristic energy scales ( "breaks" ) 
at which one loss mechanism comes to dominate over an- 
other. To fix notation, for our purposes cosmic rays are 
most usefully characterized by the propagated cosmic- 
ray flux (strictly speaking, specific intensity) (j){E) = 
dAfcr/dA dt dVl dE per unit energy interval. For all but 
the most ultra-high energies, cosmic rays are observed 
to be isotropic, in which case the flux is related to the 
cosmic-ray number density n via 4:TT(f){E) — v dn/dE. 
Here E is the cosmic-ray kinetic energy; the total rela- 
tivistic energy is thus i^tot = E+mc^. Relativistic energy 
and momentum are related by E^^^^ = [cpY + {mpC^Y , 
and v/c = cp/Etot- Using these, it follows that the flux 
per unit energy is equal to the particle number density 
per unit momentum: 47r^(i?) = dn/ dp. Hence, a number 
spectrum dn/dp that is a power law in p gives a flux with 
the same power-law of p{E). 

As a result, we characterize possible propagated proton 
spectra with a piecewise power law in relativistic momen- 
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Fig. 2. — The various cosmic-ray spectra used in this paper 
dotted curve is the leaky box propagated spectrum (710W = 0.8) 
solid curve is the broken power law spectrum with 7ic,w = —2.0 
dashed curve is the carrot spectrum with a = —4.3, / = 0.01 
above 0.2 GeV all three spectra are identical; vertical dashed lines 
represent 2 MeV and 10 MeV low energy cutoffs. 



tum p{E): 



p(Ei) 



iE2)j ' 



0, 



ifE > E2 



iiE,^t<E<E2- 



HE < K 



(8) 

Here Ei = 1 GeV is the arbitrary energy at which 
the f l ux is normalized to fit observations; following 
iMoril ()1997[ ) we take (/)p(£'i) = 0.22 cm'^ s'^ si'^ 
GeV~^. This and the observed (high-energy) spectral 
index 7iiigh ~ —2.7 fixes the high-energy region of the 
spectrum. 

The low-energy portion of the spectrum is crucial for 
this paper, and we take the high/low energy break to be 
E2 — 0.2 GeV, which is roughly where ionization losses 
begin to dominate diffusion and/or escape losses. The 
power law index for energies below E2 is 710W , and eq. ([5]) 
is arranged to guarantee that the flux is continuous across 
this break. Finally, an effective low energy cutoff, i?cut , is 
chosen, below which the flux is zero. Despite the fact that 
the flux will change as cosmic rays travel into a cloud, 
we assume a steady state such that the spectrum is the 
same everywhere. We also ne glect the p ossible effects of 
self-confinement proposed by iPadoan fc Scald (j2005i) . in 
which magnetohydrodynamics can spatially confine cos- 
mic rays to given regions due to changes in the ambient 
density. 

In the case of the propagated spectrum, the momen- 
tum index 710W = 3-1- 7sourcc — 0.8, which corresponds 
to a source spectrum with q{p) cx p~^'^ and propagation 
dominated by energy losses (the "thick-target" approx- 
imation). This spectrum is shown as the dotted curve 
in both Figures [T] & [2] The two vertical dashed lines 
in Figure [2] represent low energy cutoffs at 2 MeV and 
10 MeV. These were chosen because cosmic rays with 
these energies have ranges roughly corresponding to the 
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column densities of diffuse a nd dense clouds, respectively 
(jCravens fc Dalgarnol Il978l ). Following this reasoning, 
cosmic rays with E < 2 MeV should not penetrate dif- 
fuse clouds, and so will not contribute to the ionization 
rate there. Likewise, cosmic rays with energies below 10 
MeV will not affect the ionization rate in dense clouds. 

A nother spectrum we consid er is modeled af- 
ter iMeneguzzi et"al] ()197lD and iMeneguzzi fc ReevesI 
(|1975allbr i who added a second sharply-peaked compo- 
nent - dubbed a "carrot" - to the propagated spectrum 
to give a high flux at energies of a few MeV. The physical 
reasoning behind such a component is that in addition 
to the propagated spectrum, there is some local source 
of cosmic rays formed in weak shocks, represented by a 
steeper power law. This component is given by 



where / sets the flux of the carrot component to be some 
fraction of the propagated spectrum at E2, and the to- 
tal cosmic-ray spectrum is taken to be the sum of the 
propagated and carrot components. To ensure that the 
carrot does not conflict with observations at high ener- 
gies, / should be relatively small (< 0.1) and a must be 
less than 7high- 

In addition to the cosmic-ray spectra proposed above, 
we also consider several which have been used in the 
past for similar calculations (see §1 and Figure [T]). De- 
termining the ionization rate produced by these spectra 
allowed us to check that our numerical integration code 
was working properly, and to see exactly what energy 
range of cosmic rays is responsible for the high ioniza- 
tion rate inferred from H;^. 

5. IONIZATION ENERGETICS: THEORETICAL 
ESTIMATES 

Given a cosmic-ray spectrum and relevant ionization 
cross section, the cosmic-ray ionization rate is readily 
calculable. Namely, the ionization rate of species X due 
to cosmic rays is given by: 



Cx = 4^ex(i + qro]) 



0p(S)fTion(S) dE, (10) 



where (t>p{E) is the flux of cosmic-ray protons as a func- 
tion of kinetic energy, aioniE) is the ionization cross sec- 
tion of atomic hydrogen as a function of kinetic energy, 
Gjj^ is a coefficient accounting for ionization by heav- 
ier cosmic ray nuclei, and converts between the pri- 
mary ionization rate of atomic hydrogen (Cp), computed 
by the integral, and the total ionization rate for a given 
species X. This conversion factor includes ionization due 
to secondary electrons produced in the initial event, and 
accounts for the difference in the ionization cross section 
between H and X (Cx = ^xCp)- The coefficients for atomic 
(^h) and molecular (£2) hydrogen a re 1.5 and 2.3, re- 
spectively (jGlassgold fc Langerlll974D . For the ionization 
rate calculation GJj^] = 0-5 (The coefficient Gn changes 
based on context, and is labeled with a subscript indi- 
cating which equation it applies to, e.g. n = 10 in this 
case . See the Ap pendix for a more detailed discussion). 
The iBethi (|1933D cross section for primary ionization of 
atomic hydrogen 



CTion=27r(0.285)- 



2Ry /32 



In 



0.048(1 - /32)Ry 



TABLE 1 

Cosmic- Ray Ionization Rates (lO"-"^^ s~^) 





C2 


C2 




Ecut = 2 MeV 


Ecut = 10 MeV 


Spectrum 


(diffuse) 


(dense) 


Propagated^ 


1.4 


1.3 


Broken Power Law^ 


36 


8.6 


Carrot" 


37 


2.6 


Hayakawa ct al. (1961) 


165 


96 


Spitzer & Tomasko (1968) 


0.7 


0.7 


Nath & Bicrmaim (1994) 


260 


34 


Kneller et al. (2003) 


1.3 


1.0 


Ip & Axford (1985)'' 


3.6 


2.7 


Herbst & CuDoen (2006) 


0.9 


0.9 


Observational Inferences 


^ 40^^ 


~3d 


Note. — Ionization rates 


calculated are 


for molecular hy- 



drogen due to a spectrum of cosmic-ray protons and heav- 
ier nuclei with abundances greater than 10 with respect 
to hydrogen. Factors such as the 5/3 and 1.89 used by 
[Spitzor & Tomasko (1968) have been removed to calculate the 
primary ionization rate due to protons, which is then multi- 
plie d by 1.5 to account for th e heavy nuclei (1 + QToI), and 
2.3 I IGlasseold fc LaneerllTgT^ to find the H2 ionization rate. 

tu dy'' via fitting function of 

pOO^ )'^ llndriolo et all pOOTi ): 



this 
IValle et aU 
Ht' TraiTder Tak 



= 1.23x10-^"— (6.2 + log. 



van Dishoecg l l2000t) : Hl3cO+ 



3lO 



1-/32 



0.43/32 cm^ (12) 



(11) 



is used, where /3 = ?;/c is the velocity of the cosmic ray in 
units of the speed of light, Ry = 13.6 eV is the hydrogen 
binding energy, and Z is the cosmic-ray charge. Because 
the ionization cross section is well determined, variation 
of the cosmic-ray spectrum must be used to match the 
ionization rate inferred from observations. 

Performing a numerical integration'' of eq. (jlOp using 
the cross section from eq. ((T^ and the propagated spec- 
trum (710W = 0.8) form of eq. ^ with Scut = 2 MeV pro- 
duces a cosmic-ray ionization rate of C2 = 1.4x 10~'^ s~^, 
about 30 times smaller than the value inferred from H;^ . 
This large discrepancy shows that a simple cosmic-ray 
spectrum based on the propagation of the source spec- 
trum resulting from strong shocks is inconsistent with the 
ionization rate inferred from H;^ . To reproduce the obser- 
vational results then, we maintained the well-constrained 
high energy behavior of the cosmic-ray spectrum, but 
varied 710W despite the fact that this removes the ini- 
tial physical motivation for the low-energy portion of 
the spectrum. After several trials, we found that with 
Ecut = 2 MeV and 710W = —2.0 (shown as the solid curve 
in Figure [2]), the above calculation gives an ionization 
rate of C2 — 3.6 x 10^'^ s^'. However, when ii^cut = 10 
MeV is used to account for dense clouds, the calculated 
ionization rate is C2 = 8.6 x 10~^^ s~\ a bit larger than 
inferred values f WUliams et aL 1 9981 : iMcCafl et al.iri999l : 
Ivan der Tak fc van Dishoecki200T]f) . 

We next attempted to reproduce the inferred ion- 
ization rate by using several cosmic-ray spectra in 
the literature. These include spectra prev iously used 
to c a lculate light eleme nt abundances (jValle et al.l 
I2OO2I : iKneller et al.l l2003f l. desorption from interstel- 

Integration was performed using the qromb, trapzd, and poli nt 
routines of Numerical Recipes in FORTRAN llPress et al.ll 19921 ). 
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lar ice s (iHerbst fc Cupped [2006D. and the ioni zation 
rate^ (iHavakawa et all Il96ll : ISpitzer fc Tomaskol 119681 
iNath fc Biemiannlll994) . The results of our calculations 
using these spectra are shown in Table [1] along with the 
results from the 3 spectra proposed in this paper. It is 
interesting that none of the previous spectra are capa- 
ble of reproducing the ionization rate in diffuse clouds to 
within even the correct order of magnitude, thus high- 
lighting the need for the present study. 

To match the ionization rate in both diffuse and dense 
clouds, we then turned to the carrot spectrum, which, 
as mentioned above, must rise faster than (f)p oc p^^'^ to 
low energies. Choosing / = 0.01 and a = —4.3 (these 
values optimize the ionization and light element results 
as discussed in §6.1), and using ii^cut = 2 MeV gener- 
ates an ionization rate of 3.7 x 10^^^ s^^. The carrot 
spectrum with these parameters is shown in Figure [2] as 
the dashed curve. Changing the low energy cutoff to 
10 MeV to simulate a dense cloud environment predicts 
^2 = 2.6 X 10""'^^ s~^, also in accord with observations. 
This demonstrates that the steeper slope of the carrot 
component is better able to reproduce the roughly 1 or- 
der of magnitude difference in the ionization rate between 
diffuse and dense clouds. It is also interesting to note 
that the large majority of ionizing cosmic rays have ki- 
netic energies between 2 and 10 MeV: ~ 95% in the case 
of the carrot spectrum and ^ 80% for the broken power 
law. 

6. OTHER OBSERVABLE SIGNATURES OF LOW-ENERGY 
COSMIC-RAY INTERACTIONS 

As stated in §1, cosmic rays produce light elements and 
gamma rays via spallation and the excitation of nuclear 
states, respectively. Like the ionization rate, the pro- 
duction rates of these processes can be computed using 
eq. (fTO|) . In these cases, (Tion is replaced with the relevant 
cross section for each process, = 1, and G [jq] = 0. 

6.1. Light Elements 

To calculate the total production rates of the light el- 
ement species ^Li, ^Li, ^Be, ^°B, and ^^B (often collec- 
tively referred to as LiBeB), 32 reactions were considered. 
These include the spallation (fragmentation) reactions 

[p, a] -I- [C, N, O] [LiBeB] H , which make aU of the 

LiBeB nuclides, as well as the fusion reaction a + a 
^'^Li -|- . . ., which can only make the lit hium isotopes. 
Tabul ated cross sections were taken from iRead fc Violal 
([l98l . and for energies above ~ 50 MeV nucleon 
the a -\- a pr o cesses were supplemented with data from 
iMercer et al.l ()2001[ ). In the case of all a particle pro- 
cesses, the fluxes of the cosmic-ray spectra were reduced 
to 9.7% of the fluxes used in the ionization calculations 
because of the relative solar abundance b etween helium 
and hydrogen (|Anders fc Grevessi Il989f) . Using these 
cross sections and the spectra from §5, we calculated 
the present-day instantaneous production rate of each 
species from each process. 

Of course, the observed light element abundances 
are the result of cosmic-ray interactions with the ISM 
throughout the history of the Galaxy, meaning that they 

^ In these cases we use the same coefficients and low-energy 
cutoffs as for our proposed spectra, so our results differ slightly 
from those of the original respective papers. 



are dependent on the cosmic-ray history and chemical 
evolution of the Galaxy. A full calculation of these ef- 
fects and a comparison with LiBeB abundance evolution 
as traced by Galactic stars is a worthy subject of fu- 
ture work, but is beyond the scope of this paper. To 
estimate the accumulate d LiBeB abundance s, we follow 
the original approach of iReeves et al.l (|1970f ) to roughly 
quantify the solar LiBeB abundances expected from spal- 
lation processes with our trial spectra. Our estimate 
assumes that both the cosmic-ray spectrum and CNO 
abundances have remained constant throughout the his- 
tory of the Galaxy. Also, we assume that once created, 
the light element isotopes are not destroyed. As we know 
that light elements are destroyed ("astrated") in stellar 
interiors, this assumption leads to further uncertainty 
in our model. In addition to cosmic-ray spallation, ''Li 
and ^^B are produced by other mechanisms: neutrino 
spallation processes in Type II SNe for both isotopes 
(iDearborn et al.lll989l : IWooslev et anil990( ). and also pri- 
mordial nucleosynthesis in the case of ^Li. Both of these 
mechanisms will contribute to the observed abundances, 
but we have chosen to omit their effects with the un- 
derstanding that our ''Li and ^^B abundances should be 
lower than the net Galactic levels. These effects both add 
to and subtract from our estimate based on a constant 
production rate. Based on mor e detailed models which 
include these effects ([Fi elds & Olivd 119991: iFields et al.l 
I2OOOI : iProdanovic fc Fields. ,2006} ) we expect our calcu- 
lations of the absolute abundances to be accurate only 
to within factors of 2-3. Our results appear in Table 
[21 along with solar-system light element abundances as 
measured from meteorites and the solar photosphere. 

As seen in Table [21 the conventional propagated spec- 
trum reproduces each of the ^Li, ^Be, and ^°B abun- 
dances and their ratios well, to within 10-30%, while 
severely underproducing ^Li and ^^B. This well-known 
pattern is chara cteristic of cosmic-ray nucleosynthesis 
predi ctions (e.g.. lFields et al.ll2001l : IVangioni-Flam et al.l 
I2OOOI and references therein) and follows our expecta- 
tions given the omission of non-cosmic-ray ^Li and ^^B 
processes mentioned above. However, as we have shown, 
this spectrum leads to an ionization rate which falls far 
short of that required by data. 

Turning to the spectra with low-energy enhancements 
and associated high ionization rates, the LiBeB produc- 
tion presents a mixed picture. Here, we focus on the 
species with only cosmic-ray sources: ^Li, ^Be, and ^"B. 
In Table [2] we see that the carrot spectrum reproduces 
^°B quite well, and overproduces ^Li and ^Be by factors 
of 1.7 and 1.3, respectively, still well within the uncer- 
tainties of our rough calculation. For the broken power 
law spectrum, however, all of these absolute abundances 
exceed observations by factors of about 2-5. 

It is quite striking that there is only a factor of a 
few difference between the LiBeB abundances predicted 
by the propagated and carrot/broken power law spec- 
tra compared to the factor of about 30 difference for 
the ionization rate. This is due to two properties of the 
LiBeB production cross sections. First, most of the cross 
sections have low-energy thresholds in the tens of MeV, 
meaning that the high flux in the few MeV range has no 
effect on LiBeB production. Second, the cross sections 
do not fall off steeply as energy increases, so the cosmic- 
ray flux in the hundreds of MeV range (where all three 
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TABLE 2 
Light Element Abundances 



Ratio 



Solar System^ I Carrot Broken Power Law Propagated 



10"' X ^Li/H 


1.5 


2.5 


8.2 


1.3 


lOio X ■^Li/H 


19 


5.8 


18 


1.9 


10"' X 9Be/H 


0.26 


0.35 


0.59 


0.33 


10"' X "'B/H 


1.5 


1.4 


2.5 


1.3 


lO^" X "B/H 


6.1 


3.2 


6.4 


2.8 


e Li/9 Be 


5.8 


7.1 


13.9 


4.0 


iOB/9Be 


5.8 


4.0 


4.3 


3.9 



Note. — For all three spectra calculations were done using E^ut = 2 MeV. 
Calculated values were found by integrating the instantaneous rates over 
10 Gyr. 

^ abundances from I Anders &: Grevess^ II1989I) 



spectra are identical) is important. This contrasts with 
the case of ionization where cosmic rays with the lowest 
energies above the cutoff dominate. 

However, much more significant than the absolute 
abundances are the ratios of these isotopes, which ef- 
fectively remove the systematic uncertainties in the ab- 
solute abundances introduced by the simplicity of our 
model and directly reflect the shape of the cosmic-ray 
spectrum. While both of the low-energy-enhanced spec- 
tra underestimate the ^°B/^Be ratio by about a factor 
of 1.5, the carrot spectrum does a much better job of 
reproducing the ^Li/^Be ratio; it overestimates the ratio 
by a factor of only 1.2 compared to the broken power 
law's 2.4. 

This success of the carrot spectrum is not surprising 
though, as we chose the input parameters to best re- 
produce the observed ionization rates and light element 
ratios. These "optimal" parameters were found by using 
various combinations of / and a to compute the ioniza- 
tion rate with a 2 and 10 MeV cutoff, and the ^Li/^Be 
and ^°B/^Be ratios. Figure [3] is a plot in (/, a) space 
where the contours represent deviations of 10% and 25% 
from inferred values of ^2 in diffuse and dense clouds, 
and from measured meteoritic LiBeB ratios. It is clear 
from Figure [3] that there is an overlapping region around 
/ ^ 0.01 and a ^ —4.5 where C2 for diffuse and dense 
clouds and ^Li/^Be are all within 25% of observed val- 
ues. In making the diffuse cloud ionization rate as close 
to 4 X 10~^^ s~^ as possible, we chose / = 0.01 and 
a — —4.3 (indicated by the triangle in Figure [3]) as the 
parameters for the carrot spectrum. 

While C2 and ^Li/^Be can be matched well, no combi- 
nation of / and a is able to successfully reproduce the 
^°B/^Be ratio to within 25% simultaneously with any of 
the other observables. Almost the entire range of Fig- 
ure [3] is within the 50% contour of ^^B/^Be though, so 
our carrot spectrum is not completely out of the ques- 
tion. Indeed, the carrot spectrum predicts almost the 
same ^°B/^Be ratio as does the propagated spectrum, so 
the introduction of the carrot leaves the agreement with 
solar system data no worse than in the standard case. 

6.2. Gamma Rays 

In computing the production rates of gamma-rays, 6 
total reactions were used. These include [p, a] -I- ^^C 
12c* ^ 12c + 74.44; [P, a] + ^ 160* ^ + 76.13; 
and [p, a] + i^O '^'^C* + a ^"^C + a + 74.44, where 
the de-excitations of i^C* and ^^O* produce 4.44 MeV 
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Fig. 3. — Plot in (/, «) space showing how well various combi- 
nations reproduce the observed ionization rates and light element 
ratios. Thicker lines represent 10% deviation from observed val- 
ues, and thinner lines show 25% deviation. Different styles have 
the following meanings: solid - 1^2 in diffuse clouds; dashed - 1J2 in 
dense clouds; dotted - ®Li/®Be; dash-dot - '^"B/^Be. The triangle 
shows the parameters chosen for the carrot spectrum; / = 0.01, 
a = -4.3 

and 6.13 MeV gamma rays, respectively. Cross sections 
for all of the above processes come from iRamatv et al.l 
(fl973) (and references therein). Along a line of sight 
with hydrogen column density iVjj, the gamma-ray line 
intensity is 



iVn / HE)<J^iE) dE, 



(13) 



with units [photons cm^^ sr^^], and where for each 
line the product (f)a represents a sum over appropri- 
ately weighted projectiles and targets. If all diffuse 
clouds experience the same cosmic-ray flux and the col- 
umn of such clouds through the Galactic plane is about 
iVn ~ 10^^ cm^^, our calculations predict that the Galac- 
tic plane should have the diffuse 7-ray fluxes shown in 
Table [3l Assuming the intensity in eq. (|13p is uniform 
in the Galactic plane, we integrated over the solid an- 
gle within 1^1 < 30°, \b\ < 10° to find the total fiux in 
the central radian of the Galaxy (a typical region over 
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which diffuse 7-ray fluxes are quoted, with units [cm ^ 
s~^ rad^-'^]). Also shown in Table [3] is the large-scale 
sensitivity of the INTEGRAL spectrometer at ~ 5 MeV 
()Teegarden fc Watanabel[2b06[ ). Our predicted fluxes are 
below currently available detector limits of INTEGRAL. 
Thus the presence of low-energy cosmic-rays sufficient to 
give the ionization levels required by H;^ does not violate 
gamma-ray constraints. 

Indeed, the gamma-ray line predictions in Table [3] lie 
tantalizingly close to present limits. While this does 
not provide a test of the predicted cosmic-ray spectra 
at present, it may be possible that INTEGRAL itself, 
and certainly the next generation gamma-ray observa- 
tory, will have the ability to detect these lines. In any 
case, the results show that our proposed spectra are not 
inconsistent with observations. 

6.3. Energetics 

Similar to the calculations in §3, we can determine the 
energy budget for a given cosmic-ray spectrum. Unlike 
the previous calculation though, in this case the shape 
of the spectrum is important as we compute the energy 
loss rate of all of the cosmic rays in the spectrum. This 
is done via the usual Bcthc-Bloch expression for energy 
loss dE per unit mass column dR — pdx — pvdt: 



dE _ ^TTzZ'^e^ 
dR A{m)mev'^ 



In 



f3' 



(14) 



which is closely related to the ionization cross section 
above (eq. [12]) . Here we u se z = Z = A = 1 ( m) = 
1.4TOp, and / = 13.6 eV (see iProdanovic fc Fieldd (pOOj ) 
for a complete description of the variables involved). The 
rate of cosmic-ray energy loss per unit mass of neutral 
hydrogen is 



L 



CR 



ME) 



Bio 



dE 
dR 



dE, 



(15) 



where G|i5] =0.1 (see the Appendix). Again using a 
Galactic hydrogen mass of Mh = 5 x IO^Mq, the carrot 
and broken power law spectra require energy inputs of 
LcR = 0.18 X 10^^ erg (100 yr)-i and 0.17 x lO^^ erg 
(100 yr)~^, respectively. These represent large fractions 
of the total mechanical energy released in SNe. Like in 
§3, they are als o consistent with th e 0.5 x 10^^ erg (100 
yr)~^ found by [Fields et al.l (|2001f) which accounted for 
both energy needed for ionization of the ISM and escape 
from the Galaxy. Finally, we note that the large amounts 
of energy and high acceleration efhciencies required may 
he lp to resolve the super bubble "energy crisis" described 
bv lButt fc B"^^ (pOOSh . 

Beyond the total input energy requirement, each 
cosmic-ray spectrum will also have a particular energy 
density and pressure. Energy density can be calculated 
from 



= 4.(1 + qm) ^^^^^ 



(16) 



where v{E) is the velocity and C ^YE\ ^ '^•^-^ (^®^ ^P" 
pendix), and pressure from 



CR = + qni) ME)p{E)dE, (17) 



■E 0.01 




0.1 1 
Energy (GeV/nucleon) 

Fig. 4. — Contribution to the energy density of cosmic rays 
as a function of kinetic energy per nucleon. As in Figure [2] 
the dotted curve is the leaky box propagated spectrum, the sohd 
curve is the broken power law spectrum, the dashed curve is the 
carrot spectrum, and the vertical dashed line shows the 2 MeV 
low energy cutoff. The vertical axis is given by Ede/dE, where 
de/dE = A-KE<f)p{E)/v{E). 

where p{E) is again relativistic momentum, and Gjjy] = 
0.42. Performing these calculations with Ecut = 
2 MeV results in energy densities of 0.77 eV cm~^ and 
0.89 eV cm-3^ and pressures (Pcr/^b) of 4300 K cm'^ 
and 5200 K cm~'^ for the carrot and broken power law 
spectra, respectively. Both pressures are in rough accord 
with the average thermal pressure in t he diffuse ISM of 
P/fcg ) = 2700 K cm"3 reported by IJenkins fc Tripd 
2OOI . The energy densities in both spectra, however, 
are ab out one half of the value reported by IWebbeH 
(1998*). While this result may at first seem counterin- 
tuitive, it is best clarified graphically by Figure H) Here, 
it is shown that cosmic rays with energies between about 
0.1 GeV and 10 GeV completely dominate in contribut- 
ing to the energy density. This was pr eviously demon - 
strated by the analogous plot (Fig. 7) in Webber (1998), 
from which the author concluded that low energy compo- 
nents, such as those proposed in this paper, would have 
little effect on the cosmic-ray energy density. As for why 
our spectra have lower energy densities, this is almost 
entirely dependent on the flux normalization at higher 
energies. At about 1 GeV the fluxes in our spectra are 
ab out one half tha t of the local interstellar spectrum used 
by IWebbej ()1998D , thus resulting in the corresponding 
factor of 2 difference in energy densities. 

6.4. Cloud Heating 

One further effect that cosmic rays have is to heat 
the ISM via energy lost during the ionization process. 
On average, 30 eV are lost by a cosmic ray during 
each ionization event ([Cravens k Dalgarnd[1978[) . Using 
^2 = 4 X 10~^^ s~^ and the corresponding = 1.5^2/2.3 
and assuming that the number density of atomic hydro- 
gen is roughly equal to that of molecular hydrogen, and 
that all of the lost energy eventually goes into heating, we 
find a heating rate due to cosmic rays of Fcr = 3 x 10~^^ 
erg s~^ (H atom)^^. This can be compared to the heat- 
ing rate due to the photoelectric effect calculated by 
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TABLE 3 

Diffuse Gamma- Ray Flux from the Central Radian'' (10" 

H,AD~^) 



Energy 


INTEGRAL sensitivity 


Carrot 


Broken Power Law 


Propagated 


4.44 MeV 


10 


1 3.0 


8.3 


0.9 


6.13 MeV 


10 


1 2.4 


5.9 


0.4 



Note. — Predicted fluxes for the 4.44 and 6.13 MeV 7-ray lines using our carrot, 
broken power law, and propagated spectra. All calculations were done using -Bcut = 
2 MeV. Also shown are the most directly comparable sensitivities of the INTEGRAL 
spectrometer given by Teegardcn & Watanabc ( 2006) 

^ For the central radian wo considered uniform emission within \l\ < 30° and |6| < 
10° 



iBakes fc Tiele^ (Il99l of Tpe = 1.5 x IQ-^^ erg s'^ 
(H atom)~^ for the diffuse cloud sight line toward Q 
Oph. The heating rate caused by cosmic-ray ionization 
is about 5 times smaller than the heating rate due to the 
photoelectric effect, demonstrating that even with such 
a high ionization rate, cosmic rays do not significantly 
alter the heating rate in diffuse clouds. This large dif- 
feren ce in heating r ates i s further illustrated by Figure 
10 of lWolfire etall (|2003l ). Because the high flux of low 
energy cosmic rays we use will not dominate cloud heat- 
ing, our spectra do not imply cloud temperatures that 
are inconsistent with observations. 

7. DISCUSSION 

7.1. Cosmic-Ray Spectra 

As described in §4, the spectrum with seemingly the 
best physical motivation is one that arises from the prop- 
agation of particles accelerated by strong shocks in su- 
pernova remnants. This spectrum follows a p~^ '^ rela- 
tionship above a few hundred MeV, matching observa- 
tions, and a relationship below a few hundred MeV 
in the ionization-dominated regime. A similar behavior is 
apparent in the spectra u sed byiHavakawa et al. (1961), 
iSpitzer fc Tomaskol (|1968D . and lHerbst fc Cuppen (.2006i) 
(see Figure [T]), except they follow power laws closer to 
at low energies. Even the more sophisticated models such 
as those considering re-acceleration (Shibata ct al. 2006) 
or distributed cosmic-ray s ources and a Galactic wind 
(jLerche fc Schhckeiseijll982l ) follow these general power 
laws. Because they d ecrease at low eiiergy th ough, all of 
these spectra (except iHavakawa et al.l (|1961t ) which does 
not begin decreasing until £^ < 10 MeV) are unable to 
provide enough flux at the energies where ionization is 
most efficient, and thus cannot generate the high ion- 
ization rate inferred from H^. It seems then that the 
propagation of cosmic rays accelerated by SNR's (with 
test-particle power-law spectra) cannot explain the high 
flux of low energy cosmic rays necessary in our spectra. 
However, plentiful evidence supports the idea that high- 
energy (> 1 GeV) Galactic cosmic rays do originate in 
supernovae: synchrotron emission in supernova remnants 
indicates electron acceleration and strongly suggests ion 
acceleration, and Galaxy-wide cosmic-ray energetics are 
in line with supernova expectations and difficult to sat- 
isfy otherwise. For this reason we retained the propa- 
gated cosmic-ray spectrum to represent the supernova- 
accelerated Galactic cosmic-ray component at high ener- 
gies (> 1 GcV), but also considered additional cosmic- 
ray components which dominate at low energies where 
the ionization efficiency is high. 



While our carrot and broken power law spectra do not 
follow the conventional strong-shock thick-target (j) oc 
p"'^ relationship in the low-energy regime, recent studies 
are beginning to find possible sources for the proposed 
high flux of low energy cosmic rays. Weak shocks will 
accelerate cosmic rays with preferentially steeper power 
laws (see §4), and may be ubiquitous in the ISM, caused 
by star forming regions, OB associa- t ions, and even low 
mass stars like the sun. IStone et al.l (|2005| ) investigated 
the flux of termination shock protons at the heliosheath 
using the Voyager spacecraft. Their study showed that 
at low energies, 0.5 MeV < -B < 3.5 MeV, (jjp oc E'^-'^ 
(oc p~^'^), and that from about 3.5-10 MeV, this relation- 
ship steepened to cjip oc E~'^ (oc p~^)- These so-called 
"anomalous" cosmic rays (ACRs) clearly demonstrate 
that real sources exist in nature which can produce a high 
particle flux at low energies, though these measurements 
and the anomalous cosmic rays themselv es are located 
in the heliosphere, not interstellar space. IScherer et al.l 
()2008| ) modeled the effects from the astropauses of all F, 
G, and K type stars in the Galaxy to find a power of 
2.2 X 10"*^ erg (100 yr)-i in ACRs. This accounts for 
only about 10% of the power needed to produce the ion- 
ization rate inferred from H^. However, their analysis 
did not include the effects of winds from the much more 
luminous O and B stars. To our knowledge, no study 
has computed the interstellar cosmic-ray spectrum aris- 
ing from the ACRs of all stellar winds in the Galaxy, but 
it may indeed be an important contribution to the flux 
of low energy cosmic rays. 

Another intriguing possibility is that supernova 
shocks are considerably more efficient at low-energy 
particle acceleration than simple tcst-particle results 
would indicate. Indeed, theoret ical non-linear shock 
acceleration calculations (e.g. , iKang fc Joned Il995t 
iBerezhko fc EllisonI Il999t iBlasil 1200^ do predict that 
low-energy particles have higher fluxes than in the test- 
particle limit. This result goes in the right direction 
qualitatively. However, published source spectra we are 
aware of do not appear steep enough at low energies - 
particularly after ionization losses are taken into account 
in propagation - to reproduce the observed ionization 
rate. It remains an interesting question whether more de- 
tailed nonlinear calculations, focussing on the low-energy 
regime, might provide a solution to the ionization prob- 
lem; in this case, ionization rates and gamma-ray lines 
would become new probes of feedback processes in su- 
pernova remnants. 

In addition to protons and heavy nuclei, it has been 
proposed that cosmic-ray electrons may make a signifi- 
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cant contribution to the ionization rate. IWebbeil ()1998D 
showed that the local interstellar spectrum of cosmic- 
ray electrons produces a primary ionization rate of ~ 
2 X 10~^^ s~^ when considering energies above ~ 2 MeV. 
While this is roughly equal to the ionization rate inferred 
for diffuse clouds at that time, electrons only account 
for about 10% of the primar y ionization rate inferred 
from (jlndriolo et al.|[2"007( ). As a result, we have ne- 
glected the effect of cosmic-ray electrons in the present 
study. It is worth noting, however, that low-energy 
cosmic-ray electrons are probed by very low-frequency 
radio emission, and indeed cosmic-ray electron emission 
provides a major foreground for present and future fa- 
cilities aimed at the measurement of cosmological 21-cm 
emission from high-redshift sources, including L OFAR 
(iRottg cring 2003) and the Square Kilometer Array (|Beckl 
|2005| ). Such observations should provide a detailed pic- 
ture of low-energy cosmic-ray electrons, whose behavior 
can in turn be compared to the proton and nuclear com- 
ponents probed by the other observables considered in 
this paper. 

Finally, we have found that the carrot spectrum pro- 
duces ^Li/^Be in good agreement with solar system data, 
and a -"^^B/^Be ratio which is almost identical to the stan- 
dard propagated result but which is somewhat low with 
respect to the solar system value. To the extent that 
the isotope ratios are not in perfect agreement with so- 
lar system data, one possible explanation could be time 
variations of the cosmic-ray spectral shape over Galac- 
tic history. If supernovae are the agents of cosmic-ray 
acceleration, then presumably strong shocks will always 
lead to high-energy source spectra with 7sourcc ~ —2.2 as 
we have today. However, the low-energy component re- 
sponsible for the carrot derives from weaker shocks which 
in turn may reflect time-dependent properties of, e.g., 
star-forming regions. Moreover, cosmic-ray propagation 
is much more sensitive to the details of the interstellar en- 
viron ment, particularly the nature of Galactic magnetic 
fields. iPrantzos et al.l (|1993[ ) suggested that such varia- 
tions in the early Galaxy might explain the B/Be ratios 
in primitive (Population II) halo stars. Similarly, if such 
variations were present in the later phases of Galactic 
evolution (e.g., during major merger events) then it is 
possible that the propagated cosmic-ray spectra could 
have differed substantially. The resulting LiBeB con- 
tributions could alter the ratios from the simple time- 
independent estimates we have made. 

Another possible explanation to bring the theoreti- 
cal LiBeB ratios into better agreement with observa- 
tions would arise if the LiBeB isotopes suffer signifi- 
cantly different amounts of destruction (astration) in 
stellar interiors. Because the binding energies are in 
the hierarchy BC^U) < B{^Bc) < B(^^B), there should 
be a similar ranking of the fraction of the initial stel- 
lar abundance of these isotopes which survives to be 
re-ejected at a star's death. If the amount of ^Li de- 
stroyed is greater than ^Be, which is in turn greater than 
^"^B, then the results from our carrot spectrum may be 
correct before accounting for astration. Assuming this 
preferential destruction decreases our calculated ^Li/^Be 
and increases ^'^B/^Be, moving both closer to the mea- 
sured ratios. That s aid, conventional stellar models (e.g., 
ISackmann fc Boothrovd 1999.) a nd their implement ation 
in Galactic chemical evolution (|Alibes et al.l l200"^ find 



different, but still small, survival fractions for all iso- 
topes, < 10% for ^°B. As a result, this scenario would 
seem to require large upward revisions to the survival of 
^Be and ^°B in stars. 

7.2. Astrochemistry 

Gas phase chemistry in the ISM is driven by ion- 
molecule reactions. Photons with E > 13.6 eV are 
severely attenuated in diffuse and dense clouds, meaning 
that cosmic rays are the primary ionization mechanism 
in such environments. As a result, the cosmic-ray ioniza- 
tion rate has a large effect on the chemical complexity of 
the cold neutral medium. In fact, it has a direct impact 
on the abundances of H+, OH, HD, HCO+, and H3O+, 
to name a few molecules. This makes the cosmic-ray 
ionization rate an important input parameter for astro- 
chemical models which predict the abundances of various 
atomic and molecular species. 

However, based on the current theoretical study it 
seems that instead of having a uniform Galactic value, 
the cosmic-ray ionization rate should vary between sight 
lines. This has to do with the source behind cosmic- 
ray acceleration. While it has typically been assumed 
that low energy cosmic rays are accelerated in supernovae 
remnants, the spectra making this assumption were un- 
able to reproduce the ionization rate inferred from . 
Instead, a low energy carrot was required, most likely 
produced by particles accelerated in weaker, more lo- 
calized shocks. Unlike the SNe cosmic rays which are 
assumed to diffuse throughout the Galaxy, cosmic rays 
accelerated in weak local shocks could lead to significant 
enhancements in the local ionization rate. Observations 
of H;^ have shown that the cosmic-ray ionization rate is 
in fact variable between different diffuse cloud sight lines. 
The H2 ioniza tion rates towa.rd C P er and X Per are about 
7 X 10"^^ s-i findriolo et al.ll2007D . while 3cr upper Hmits 
toward C Oph and o Sco are as low as 1.6 x 10~^^ s~^ 
and 1.2 x 10~^^ s~\ res pectively^. Like the results of 
Ivan der Tak et al] (|2006D for dense clouds, this demon- 
strates that the cosmic-ray ionization rate can vary sig- 
nificantly between diffuse clouds as well. Also, it suggests 
that instead of searching for or adopting a "canonical" 
ionization rate, sight lines must be evaluated on a case- 
by-case basis. 

In order to test the theory that low energy cosmic rays 
are primarily accelerated by localized shocks, we propose 
two complementary observational surveys. First, the ion- 
ization rate should be inferred along several diffuse cloud 
sight lines which are surrounded by different environ- 
ments. Observations of Hj^ in sight lines near OB associ- 
ations and sight lines near low mass stars should provide 
data in regions near and far from energetic sources. We 
expect the sight lines near more energetic regions to show 
higher ionization rates than those near less energetic re- 
gions. If observations confirm these predictions, then we 
will be able to more confidently conclude that most of 
the low energy ionizing cosmic rays are accelerated in 
localized shocks. 

The second set of observations we propose examines 
the ionization rate in regions of varying density. FoUow- 

These li mits are based on ob servations performed after the 
publication of llndriolo et all 1)20071 ) . and will be described in more 
detail in a future publication 
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ing the reasoning of §4 where we assume that the lower 
energy cosmic rays do not penetrate denser clouds, the 
ionization rate should be inversely related to the cloud 
density. Observing H;^ in diffuse clouds, dense clouds, 
and in sight lines with intermediate densities should pro- 
vide us with a range of ionization rates. We can then use 
the carrot spectrum with appropriate low energy cut- 
offs to reproduce the inferred ionization rates from each 
environment, thus constraining the slope of the carrot 
component. This slope will then allow us to roughly de- 
termine the strength (or rather weakness) of the shock 
necessary to produce such a steep power law, and thus 
infer the source of the shock. 

8. CONCLUSIONS 

Three theoretical low energy cosmic-ray spectra have 
been examined, all of which are consistent with di- 
rect cosmic-ray observations at high energies. We first 
adopted the standard q oc p^^-^ source spectrum result- 
ing from cosmic-ray acceleration by strong shocks in su- 
pernovae. The propagated version of this spectrum pro- 
duced an ionization rate about 30 times smaller than 
that inferred from H|^, thus demanding that additional 
low-energy cosmic-rays be responsible for the observed 
ionization. 

We thus studied the effects of ad hoc but physically 
motivated low-energy cosmic-ray components. We found 
that the carrot and broken power law spectra could be 
fashioned so as to reproduce observed results for diffuse 
clouds. Out of these two, the carrot spectrum did a much 
better job of matching the ionization rate in dense clouds. 
Unlike the broken power law, the carrot spectrum was 
also capable of matching observed light element abun- 



dances to within a factor of 2 for the three isotopes pro- 
duced solely by cosmic-ray spallation. These results are 
well within the expected uncertainties of LiBeB Galactic 
chemical evolution. 

Predictions for the gamma-ray line fluxes for both the 
carrot and broken power law were below the limits of 
current instruments, so these spectra are not inconsistent 
with data. Indeed, our estimates are close to the detec- 
tion limits of INTEGRAL, and thus motivate a search 
for the 4.44 MeV and 6.13 MeV lines (or limits to their 
intensity) in the gamma-ray sky. 

Finally, the energy necessary to accelerate all of the 
cosmic rays in these spectra is about 0.2 x 10^^ erg (100 
yr)~^. This is a substantial fraction of the mechanical 
energy released in supernovae explosions, although based 
on our results it may be necessary that much of this 
energy come from weak shocks in order to produce a 
high flux of low energy cosmic rays. 

Together, all of these calculations demonstrate that the 
proposed carrot spectrum is consistent with astrochem- 
ical and astrophysical constraints. Whether or not low- 
energy cosmic rays take precisely this spectral form, at 
the very least this example serves as a proof by construc- 
tion that one can make cosmic-ray models which contain 
low-energy enhancements required by the high ionization 
rate inferred from H;^, while not grossly violating other 
observational constraints. This motivates future work 
which looks in more detail at the impact of low-energy 
cosmic rays. 

The authors would like to thank T. Oka for insightful 
discussions. NI and BJM have been supported by NSF 
grant PHY 05-55486. 



APPENDIX 

APPENDIX: EFFECTS OF HEAVY COSMIC-RAY NUCLEI 

For the sake of clarity, discussions of the cosmic-ray spectra in the body of the paper focused only on the proton 
spectrum. Our calculations, however, included the effects of heavier nuclei cosmic rays as a coefficient in eqs. (flOl [TSl 
\TE\ & [17|) . This appendix discusses in more detail the calculations that went into determining the coefficients G„ . 

We assume that all heavy cosmic-ray nuclei have the same spectral shape as protons, but that their fluxes are shifted 
down by their respective relative abundances to hydrogen (e.g. 0hc(£') = 0.097(/)p(i?)). With this assumption, the 
contribution to the ionization rate due to heavy nuclei can be calculated from 

qm^J2^'9r, (Ai) 

i 

where Zf is the charge which comes from eq. (fT2|) . gi is the fractional abundance with respect to hydrogen, and 
the index i sums o ver aU species with solar abundances > 10"^ (^He, ^^C, ^'^N, ^^^O, ^^Ne, ^'^Mg, ^^Si, ^^S, ^"^Fe 
(|Mever et al.lll998[ )). Performing this summation results in the value of Gj^Q] = 0.5 used in the ionization calculations 
of §5. In §6.3, however, the energy loss per unit hydrogen mass (eq. llSp is controlled by the particle energy loss per 
unit mass column dE/dR oc Z'^ /A (eq. fT4|) . This changes the heavy nuclei coefficient to 

i 

where the atomic mass, Ai, is now included because of eq. I|14p . As a result, for the energy loss rate calculation 
G[j^= 0.1. Also in §6.3, the energy density (eq. [TB)) and pressure (eq. [T7|) calculations both require the coefficient 

'3m=^=Y.^-9^■ (A3) 

i 

Here, Ai is required because E and thus p{E) are both in units of per nucleon throughout the paper. In this case, 
Gjjg]= GjYT]— 0-42 for both the cosmic-ray energy density and pressure calculations. 

TTowever, if the relative abundances of Galactic cosmic rays (GCRs) measured at higher energies ()Mever et al.lllQQSl) 
are used instead of solar abundances, the above coefficients change. This is because the abundances of most heavy 
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nuclei are enhanced in GCRs when compared to solar. For the case of ionization, Gh^ becomes 1.4, making heavy 
nuclei more important than protons. Because the integral is multiplied by (1 + Qjoj) though, the overall difference in 
the ionization rate between using solar and GCR abundances is only a factor of 2^71-5=l-6. Using GCR abundances 
to calculate the energy loss rate changes Gfjglby a negligible amount, from 0.1 to 0.11. Finally, GCR abundances only 
change G[Yg] and G[Y7]from 0.42 to 0.46 forThe energy density and pressure calculations. Despite the fact that heavy 
nuclei are measured to be more abundant in Galactic cosmic rays than in the solar system, we have chosen to use solar 
abundances in our calculations. This is because the high energy cosmic rays observed are accelerated in metal-rich 
SNRs, while the high flux of low energy cosmic rays is most likely due to weak shocks associated with low mass stars 
and the ISM. Due to this source difference, we find it justifiable to use solar abundances instead of the measured high 
energy GCR abundances. 

REFERENCES 



Alibes, A., Labay, J., & Canal, R. 2002, ApJ, 571, 326 
AMS Collaboration, et al. 2002, Phys. Rep., 366, 331 
Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Arnett, W. D. 1987, ApJ, 319, 136 

Bakes, E. L. O., & Tielens, A. G. G. M. 1994, ApJ, 427, 822 
Beck, R. 2005, Astronomische Nachrichten, 326, 608 
Bell, A. R. 1978, MNRAS, 182, 147 

Bethe, H. 1933, Handbuch der Physik, (Berlin: J. Springer), 24, 
Ft. 1, 491 

Black, J. H., & Dalgarno, A. 1977, ApJS, 34, 405 

Black, J. H., Hartquist, T. W., & Dalgarno, A. 1978, ApJ, 224, 448 

Blandford, R. D., & Ostriker, J. F. 1978, ApJ, 221, L29 

Berezhko, E. G., & Ellison, D. C. 1999, ApJ, 526, 385 

Blasi, F. 2002, Astroparticle Physics, 16, 429 

Brinks, E. 1990 in Proc. 2nd Teton Conf., Vol. 161, Interstellar 
Medium in Galaxies, ed. H. A. Thronson, Jr., & J. M. ShuU 
(Dordrecht: Kluwer), 39 
Butt, Y. M., & Bykov, A. M. 2008, ApJ, 677, L21 
Cardelli, J. A., Meyer, D. M., Jura, M., & Savage, B. D. 1996, ApJ, 
467, 334 

Casse, M., Lehoucq, R., & Vangioni-Flam, E. 1995, Nature, 373, 
318 

Cravens, T. E., & Dalgarno, A. 1978, ApJ, 219, 750 
Dalgarno, A. 2006, PNAS, 103, 12269 

Dearborn, D. S. P., Schramm, D. N., Steigman, G., & Truran, J. 

1989 ApJ, 347, 455 
Dragicevich, P. M., Blair, D. G., & Burman, R. R. 1999, MNRAS, 

302, 693 

Ellison, D. C, Drury, L. O'C, & Meyer, J.-P. 1997, ApJ, 487, 197 
Federman, S. R., Weber, J., & Lambert, D. L. 1996, ApJ, 463, 181 
Fields, B. D., Casse, M., Vangioni-Flam, E., & Nomoto, K. 1996, 

ApJ, 462, 276 
Fields, B. D., & Olive, K. A. 1999, ApJ, 516, 797 
Fields, B. D., Olive, K. A., Vangioni-Flam, E., & Casse, M. 2000, 

ApJ, 540, 930 

Fields, B. D., Olive, K. A., Casse, M., & Vangioni-Flam, E. 2001, 
A&A, 370, 623 

Geballe, T. R., McCall, B. J., Hinkle, K. H., & Oka, T. 1999, ApJ, 
510, 251 

Glassgold, A. E., & Langer, W. D. 1974, ApJ, 193, 73 
Gloeckler, G., & Jokipii, J. R. 1969, Phys. Rev. Lett., 22, 1448 
Goto, M., et al. 2008, ApJ, 688, 306 

Hartquist, T. W., Black, J. H., & Dalgarno, A. 1978a, MNRAS, 
185, 643 

Hartquist, T. W., Doyle, H. T., & Dalgarno, A. 1978b, A&A, 68, 
65 

Hayakawa, S., Nishimura, S., & Takayanagi, T. 1961, PASJ, 13, 
184 

Henderson, A. P., Jackson, P. D., & Kerr, F. J. 1982, ApJ, 263, 
116 

Herbst, E., & Cuppen, H. M. 2006, PNAS, 103, 12257 
Indriolo, N., Geballe, T. R., Oka, T., & McCall, B. J. 2007, ApJ, 
671, 1736 

Inokuti, M. 1971, Reviews of Modern Physics, 43, 297 

Ip, W.-H., & Axford, W. I. 1985, A&A, 149, 7 

Jenkins, E. B., & Tripp, T. M. 2007, in ASP Conf. Ser. 365, 
Small Ionized and Neutral Structures in the Diffuse Interstellar 
Medium, ed. M. Haverkorn & W. M. Goss. (San Francisco: ASP), 
51 

Kang, H., & Jones, T. W. 1995, ApJ, 447, 944 

Kneller, J. P., PhiUips, J. R., & Walker, T. P. 2003, ApJ, 589, 217 
Krymskii, G. F. 1977, Akademiia Nauk SSSR Doklady, 234, 1306 



Le Petit, F., Roueff, E., & Herbst, E. 2004, A&A, 417, 993 
Lerche, I., & Schlickeiser, R. 1982, MNRAS, 201, 1041 
Liszt, H. 2003, A&A, 398, 621 

McCall, B. J., Geballe, T. R., Hinkle, K. H., & Oka, T. 1999, ApJ, 
522, 338 

McCall, B. J., et al. 2003, Nature, 422, 500 
McCall, B. J., et al. 2004, Phys. Rev. A, 70, 052716 
Meneguzzi, M., Audouze, J., & Reeves, H. 1971, A&A, 15, 337 
Meneguzzi, M., & Reeves, H. 1975a, A&A, 40, 91 
Meneguzzi, M., & Reeves, H. 1975b, A&A, 40, 99 
Mercer, D. J., et al. 2001, Phys. Rev. C, 63, 065805 
Meyer. J.-P., Drury, L. O'C, & EUison, D. C. 1998, Space Sci. 
Rev., 86, 179 

Misiriotis, A., Xilouris, E. M., Papamastorakis, J., Boumis, P., & 

Goudis, C. D. 2006, A&A, 459, 113 
Mori, M. 1997, ApJ, 478, 225 

Nath, B. B., & Biermann, P. L. 1994, MNRAS, 267, 447 

O'Donnell, E. J., & Watson, W. D. 1974, ApJ, 191, 89 

Oka, T., Geballe, T. R., Goto, M., Usuda, T., & McCall, B. J. 

2005, ApJ, 632, 882 
Padoan, P., & Scalo, J. 2005, ApJ, 624, L97 

Prantzos, N., Casse, M., & Vangioni-Flam, E. 1993, ApJ, 403, 630 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. 
P. 1992, Numerical Recipes in FORTRAN: The Art of Scientific 
Computing (2nd ed.; Cambridge: Cambridge Univ. Press) 
Prodanovic, T., & Fields, B. D. 2003, ApJ, 597, 48 
Prodanovic, T., & Fields, B. D. 2006, ApJ, 645, L125 
Ramaty, R., Kozlovsky, B., & Lingenfelter, R. E. 1979, ApJS, 40, 
487 

Ramaty, R., Kozlovsky, B., & Lingenfelter, R. E. 1996, ApJ, 456, 
525 

Read, S. M., & Viola, V. E. 1984, Atomic Data and Nuclear Data 
Tables, 31, 359 

Reeves, H., Fowler, W. A., & Hoyle, F. 1970, Nature, 226, 727 
Rottgering, H. 2003, New Astronomy Review, 47, 405 
Sackmann, I.- J., & Boothroyd, A. I. 1999, ApJ, 510, 217 
Scherer, K., Fichtner, H., Ferreira, S. E. S., Biisching, I., & 

Potgieter, M. S. 2008, ApJ, 680, L105 
Shaw, G., Ferland, G. J., Srianand, R., & Abel, N. P. 2006, ApJ, 

639, 941 

Shaw, G., Ferland, G. J., Srianand, R., Abel, N. P., van Hoof, P. 

A. M., & Standi, P. C. 2008, ApJ, 675, 405 
Shibata, T., Hareyama, M., Nakazawa, M., & Saito, C. 2006, ApJ, 

642, 882 

Sodroski, T. J., et al. 1994, ApJ, 428, 638 

Spitzer, L., Jr., & Tomasko, M. G. 1968, ApJ, 152, 971 

Steigman, G., & Walker, T. P. 1992, ApJ, 385, L13 

Stone, E. C, Cummings, A. C, McDonald, F. B., Heikkila, B. C, 

Lai, N., & Webber, W. R. 2005, Science, 309, 2017 
Strong, A. W., Moskalenko, I. V., & Ptuskin, V. S. 2007, Annual 

Review of Nuclear and Particle Science, 57, 285 
Tatischeff, V., & Kiener, J. 2004, New Astronomy Review, 48, 99 
Teegarden, B. J., & Watanabe, K. 2006, ApJ, 646, 965 
Valle, G., Ferrini, F., Galli, D., & Shore, S. N. 2002, ApJ, 566, 252 
van den Bergh, S., & Tammann, G. A. 1991, ARA&A, 29, 363 
van der Tak, F. F. S., Belloche, A., Schilke, P., Gtisten, R., Philipp, 

S., Comito, C, Bergman, P., & Nyman, L.-A. 2006, A&A, 454, 

L99 

van der Tak, F. F. S, & van Dishoeck, E. F. 2000, A&AL, 358, L79 
van Dishoeck, E. F., & Black, J. H. 1986, ApJS, 62, 109 
Vangioni-Flam, E., Casse, M., & Audouze, J. 2000, Phys. Rep., 
333, 365 



13 



Vangioni-Flam, E., Casse, M., Fields, B. D., & Olive, K. A. 1996, 
ApJ, 468, 199 

Wagenblast, R., & Hartquist, T. W. 1988, MNRAS, 230, 363 
Walker, T. R, Mathews, G. J., & Viola, V. E. 1985, ApJ, 299, 745 
Webber, W. R. 1998, ApJ, 506, 329 

Williams, J. P., Bergin, E. A., Caselli, P., Myers, P. C, & Plume, 

R. 1998, ApJ, 503, 689 
Wolfire, M. G., McKee, C. F., HoUenbach, D., & Tielens, A. G. G. 

M. 2003, ApJ, 587, 278 



Woosley, S. E. 1988, ApJ, 330, 218 

Woosley, S. E., Hartmann, D. H., Hoffman, R. D., & Haxton, W. 

C. 1990, ApJ, 356, 272 
Yusef-Zadeh, F., Muno, M., Wardle, M., & Lis, D. C. 2007, ApJ, 

656, 847 



